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The breakdown of Lorentz’s and CPT invariance, as described by the Extension of the Standard 
Model, gives rise to a modification of the dispersion relation of particles. Consequences of such a 
modification are reviewed in the framework of pulsar kicks induced by neutrino oscillations (active- 
sterile conversion). A peculiar feature of the modified energy-momentum relations is the occurrence 
of terms of the form <511 • p, where <511 accounts for the difference of spatial components of flavor 
depending coefficients which lead to the departure of the Lorentz symmetry, and p = p/p, being p 
the neutrino momentum. Owing to the relative orientation of p with respect to <511, the coupling 
<511 ■ p may induce the mechanism to generate the observed pulsar velocities. Topics related to the 
velocity distribution of pulsars are also discussed. 


PACS No.: 11.30.Cp, 11.30.Er, 97.60.Gb, 14.60.Pq 

The studies of a possible breakdown of the fundamental symmetries in physics represent a very active research area. 
As suggested by Kostelecky and Samuel. String/M theory provides a scenario in which a departure of the Lorentz 
invariance might manifest [1]. Recently, these investigations have been reconsidered in the context of D-branes [2], 
Loop Quantum Gravity [3-5], Non-Commutative Geometry [6], and through the spacetime variation of fundamental 
coupling constants [8]. For modern tests on Lorentz invariance, see [7] 

According to Ref. [1], Lorentz’s invariance violation (due to non trivial solution of (open) string held theory) 
follows from the observation that the vacuum solution of the theory could spontaneously violate the Lorentz and CPT 
invariance, even though such symmetries are satisfied by the underlying theory. The breakdown of these fundamental 
symmetries occurs in the Extension of the Standard Model, and only operators of mass with dimension four or less 
[9,10] are involved in order that the standard model power-counting renormalizability is preserved. In a recent work 
by Kostelecky and Mewes [11], it has been studied the general formalism for violations of Lorentz and CPT symmetry 
in the neutrino sector. The generalized equation of motion for free fermions (neutrinos in our case) is given by (we 
shall use natural units c = 1 = K) 


{iV v AB d v -M AB )^ B = 0 , ( 1 ) 

where the spinor ipg contains all the fields and their conjugates, the indices A and B range over all 2N possibility 
{fif C }) being / = e, pt, t, ... the neutrino flavors and f A = CabIb■ C is the symmetric matrix with non zero 
components C//c = 1. V U AB and Mab are 4x4 matrices in the spinor space, and can be decomposed using the basis 
of 7 matrices 

r AB = 7AB + c Ab7/x + ^Ab757m + e AB + */ab 75 + T^AB^A/i 1 (2) 

and 

M A b = rriAB + im 5A Bl5 + a l _ lA B'y' J ' + b^ABls!* 1 + ^ H AB a ^ ■ ( 3 ) 

The coefficients cl^abi b^AB,c^ B , ... are constants and in general they are flavor depending, and <7^ = A7C7"]. 
C AB’ ^AB’ an d ^ab preserve the CPT invariance, while a^AB, b^AB, e^AB, f^AB violate CPT and Lorentz invariance. 
Finally m and 7715 are Lorentz and CPT conserving. 

The time evolution of neutrinos is governed by the effective Hamiltonian 
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where (cl)^ = (c + d)^ b and (ol)^ = (a + 6)(( b , the vector (e+)" = -^(0,61 + ie 2 ) is related to the helicity of 
neutrinos (ei and e\ are two real vectors), while ( fh 2 ) a b is related to the usual neutrino masses. Details of properties 
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of the coefficients entering in the effective Hamiltonian can be found in [11]. It is worth to note that Eq. (4) involves 
only the coefficients a, c, g and H , since the other coefficients breaking the Lorentz invariance can be removed, at the 
leading order, with an appropriate redefinition of fields. 

Bounds on parameters entering in the Extension of the Standard Model which involve neutrino oscillations have been 
discussed in [11 15]. Stringent bounds on the coefficients (ol)^ and {cl) ij,v have been obtained in [11 -13] by analyzing 
various type of experiments (solar neutrino experiments [SK, SNO]; atmospheric neutrinos experiments [SK, Kam- 
LAND, LSND, K2K]; reactor experiments [CHOOZ, Palo Verde]; short-baseline experiments [CHORUS, NOMAD, 
KARMEN]; long-baseline experiments [ICARUS, MINOS, OPERA]). Estimations of the attainable sensitivities to 
the dimensionless (5cl = c£ — c£) and dimension-one ( Sa l = oJ L — a£ ) coefficients are [11] 

1(T 26 < 6dj" < 1(T 17 , 1(T 19 eV < <5a£ < 10" 9 eV. (5) 


The aim of this paper is to investigate the consequences of the Lorentz invariance breakdown in relation to pulsar 
kicks induced by neutrino oscillations in matter. The origin of the pulsar velocity is till now an open issue of the 
modern astrophysics. As follows from the observations, pulsars have a very high velocity, in comparison with the 
surrounding stars, which may varies from 450±90Km/sec up to values greater than lOOOKm/sec [16,17]. This suggests 
that nascent pulsars undergo to some kind of kick. After the supernova collapse of a massive star, neutrinos carry 
away almost all (99%) the gravitational binding energy (3 x 10 53 erg). The momentum taken by them is about 10 43 gr 
cm/sec. A fractional anisotropy of the order ~ 1% of the outgoing neutrino momenta would suffice to account for 
the neutron star recoil of 300Km/sec. Many mechanisms have been proposed to solve such a issue, but a definitive 
solution is still laking. 

An elegant and interesting mechanism to generate the pulsar velocity, which involves the neutrino oscillation physics, 
has been proposed by Kusenko and Segre [18]. Under suitable conditions, a resonant oscillation v e —» v tL T may occur 
in the region between the corresponding ( v e and z/ /tjT ) neutrinospheres. The v e are trapped by the medium (due to 
neutral and charged interactions) but neutrinos generate via oscillations can escape from the protostar being 
outside of their neutrinosphere. Thus, the surface of the resonance acts as an effective muon/tau neutrinosphere. In 
the presence of a magnetic field B, the surface of resonance is distorted with the ensuing that the energy flux turns out 
to be generated anisotropically (neutrinos generated in regions with different temperatures are emitted with different 
energies). Indeed, in a magnetized medium, the resonance condition in matter is modified by a term ~ B • p, where 
p the neutrino momentum. The relative orientation of the neutrino momentum with respect to the magnetic field 
generate the asymmetry in the neutrino emission. The case in which active-sterile neutrinos are involved has been 
studied in [19]. Papers dealing with the problem of the origin of pulsar velocities, can be found in Refs. [20-30]. 

The effective Hamiltonian (4) suggests a further mechanism to generate pulsar kicks. As we will see, terms breaking 
the Lorentz invariance give rise to a coupling of the form ~ (511 • p, where (511 is the difference of spatial components 
c 9 *, a z L referred to different flavors, and p = p/p, where p = |p|. The relative orientation of neutrino momenta with 
respect to (511 determines an asymmetry in the neutrino emission, hence can generate the pulsar kicks. This effect 
might provide a signature of a possible violation of the Lorentz invariance. 

As observed in [19], oscillations of active neutrinos could be a plausible explanation of the observed pulsar velocities 
if the resonant conversion i/ T)M <-> v e occurs between two different neutrinospheres. However, the resonant transition 
between two neutrinoshperes leads to neutrino masses which do not agree with the present limits on the masses of 
standard electroweak neutrinos. These limits do not apply to sterile neutrinos that may have only a small mixing 
angle with the ordinary neutrinos [19]. Thus we shall confine ourselves to active-sterile conversion of neutrinos. 

Eq. (4) implies that the energy of neutrinos with a given flavor is (we shall consider only the diagonal terms) 
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For our purpose, it is convenient to rewrite Eq. (7) in the form 
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and the relations p° ~ p and c l = c°* have been used. 

To estimate the anisotropy of the outgoing neutrinos, one needs to evaluate the energy flux F s emitted by nascent 
stars. According to Ref. [23], the asymmetry of the neutrino momentum is 

A^i / 0 T F^).m& 1 g 

P 3 fj F s (i9) ■ nda 9 r res ' 

The factor 1/3 takes into account of the fact that only one neutrinos species is responsible for the anisotropy (thus it 
carries out only 1/3 of the total energy). F s is the outgoing neutrino flux through the element area da of the emission 
surface. STl = H/|n| is the direction of the vector <511 = n^ / - ) — 11^, whereas n is the unity vector orthogonal to 
the element area da. The parameter g is the radial deformation of the effective surface of resonance generated by the 
term breaking the Lorentz invariance II. It shifts the resonance point r res to r(<j>) = r re s + gcos(j), with g <C r res and 
cos <j) = <5n • p. The resonance point r res is determined by the usual resonance condition (MSW effect [31]) 


2<5c2 IAj ('f‘res ) — 0? 
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Am 2 = to 2 — ?n 2 is the mass squared difference (in the notation of (13), we shall also indicate sin 29 = S 2 ). In (12), 
the potential V v . is defined as follows 

V v , = -V Pe = V 0 {3Y e - 1 + 4 Y Ve ), (14) 

K MlT = ~V^, T =V 0 (Y e -l + 2 Y Ve ), (15) 

where Y e (Y l/e ) represents the ratio between the number density of electrons (neutrinos), and 
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(16) 


m n = 938MeV is the nucleon mass and p the matter density. For sterile neutrinos one has V Vs = 0. 

The equation of evolution describing the conversion between two neutrino flavors (we consider the conversion 
v f <-> v s , f = e, p, t ) is 
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where the matrix TC is the effective Hamiltonian defined as 


n = 


v Vf - c 2 s + n ( ^ } • p s 2 <5 
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(17) 


(18) 


up to terms proportional to identity matrix 1 . 


X A comment is in order. The Lorentz invariance violation implies that the propagation of neutrinos depends, in general, on 
their flavors through the parameters Cjf, cil/i (for the moment, we shall neglect neutrino masses). The velocity is 

dE dfl 

Vf = -p- « 1 - 1- ■ 

dp dp 

Different neutrino species may have different maximum attainable velocities [32,33]. This occurs if neutrino flavor eigenstates 
are distinct from neutrino velocity eigenstates, being the two eigenstates related linearly by the mixing angle 9 V . Thus, the 
diagonal terms ■ p in (18) should be replaced by ■ pcos 2 9 V , whereas the term Tl^ l/f,l/s ^ • psin 2 9 V should appear 

in the off-diagonal. For simplicity we have taken 9 V = 0. 
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As follows from (18), the resonance condition reads 


2<5c 2 = v Vf +s n • p + <5n 0 . 


(19) 


where 


<511 • p = (i5a — 2pSc) ■ p (20) 

= |<5a — 2pSc\ cos <f>, 

and <5c = — c^ s \ 5a. = a ^ v f' > — a^ e \ A similar expression holds for (5IIo■ Notice that <5IIo(p) = <5IIo(—p) and, for 

typical values of p ~ (10 11 — 10 14 )gr/cm 3 in the protostar, (5n 0 <C 25c 2 according to constraints (5). Moreover, since 
511 • p changes sign under the transformation p —> — p, it may deform the resonance surface, a condition necessary in 
order that the asymmetry of the outgoing neutrino momenta may occur. 

Let us now evaluate p. Expanding the terms in (19) about to r — r res + gcostfi, with p —> p + 5 P , V Vf —> V Vf + 5y u 
[23], in which 

5 P = — y-^-pg cos<f> = h~ l p geos<f>, (21) 

dr v 

d InVA 

<5v„ =-:- V„ f g cos 6 = hy V Vf gcos6, (22) 

/ J v Vf J 

and using the resonance condition (12), one infers (all quantities are evaluated at the resonance) 


E 1 

V Vf hp 1 + hy * 


(23) 


where 


E = 2p\5c L \ + |<5a L |. 

Inserting Eq. (23) into Eq. (11) one gets 

Ap _ 1 E 1 

V 9 L)// r res (h p A hy ^^) 

which implies that the fractional asymmetry A p/p is ~ 1% provided 

E ~ 0.09 VyfVresihp 1 + hy 4 ). 

To compute h~ 1 + hy 1 one has to specify a model for the protostar. To this aim, we assume that the inner core of a 

protostar is consistently described by a polytropic gas of relativistic nucleon with adiabatic index T = 4/3 [34]. The 
pressure P and matter density p are related by [34,23] 

P = Kp v , (27) 

where K = T c /m n p)/ i ~ 5.6 x 10~ 5 MeV“ 4 / 3 . Here T c = 40MeV and p c = 2 x 10 14 gr/cm 3 are the temperature and 
the matter density of the core, respectively. 

The matter density p(r) can be expressed in the form [23] (see Appendix) 

p r_1 (cc) = p^^a'a; 2 + b’x + c!] , (28) 

where 

x = — , a' = (1 — p .)Ar , b' = (2 p — l)Ar , c' = 1 — p\r , (29) 

r 

and r c = 10km is the core radius. The parameter p is determined by setting p{R s ) = 0 ( R s the radius of the star) 
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Ar in (30) and (29) is given by 


. _ G n M c r - 1 n on 2 G n M c 10km 40MeV 
r—i 7 . ^ ~ 

r c p l c KV km r c T c 

where M c ~ Mq is the mass of the core (Mq is the solar mass). The temperature profile T(r) is related to matter 
density p{r) through the relation [23] (see Appendix) 


dT 2 9 kL c 
~T~ = - T P’ 

dr irr z 


(32) 


where L c is the core luminosity, L c ~ 9.5 x 10 51 erg/sec, and k = 5.6 x 10 9 cm 4 /erg 3 sec 2 ~ 6.2 x 10 56 eV 5 . Eq. 
(32) can be immediately integrate by using (28) 


T(r) = T c sj2\ c [x{x) -*(!)] + 1, 


(33) 


where 


X(x) = c' 3 x + ^ b' c' 2 x 2 + d(a! d + b' 2 ) x 3 + — (6 a'd + b' 2 ) x 4 + ( a’ d + b 12 ) x 5 H--— x 6 + x 7 , (34) 


and 


9 KL c p c p c 10km /40MeVV 

A c =--~ 1-95- T - —— 

27r T 2 r c 10 14 gr/cm 3 r c \ T c J 


(35) 


Since p ~ T (it is assumed the thermal equilibrium between neutrinos and the medium, so that the average energy 
of the emitted neutrino is proportional to the temperature at the emission point [23]) and V e ff ~ p, one can rewrite 
the inverse characteristic lengths h ~ 1 and hy 4 as h ~ 1 = hy 1 and hy 4 = h ~ 1 [23]. Eqs. (28) and (33) imply (at the 
resonance) 



d In T p{ 

^res) 

dr c p c 

dhip = _ 3 f Pc 
dr \pi.r res ) 




(2 a x res + b ) 


Xres 
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(36) 

(37) 


so that Eq. (26) reads 


E ~ 0.09 V Vf Ar x res V , (38) 

where 

77 = £ 2 ^ + 3(2p - 1) - 6(p - l)x res , (39) 

Ar 

and the parameter £ = T c /T(r res ) has been introduced. Eq. (39) gives a constraint on possible values of p, e and 
X r es • 

The Lorentz and CPT symmetry breakdown is relevant (on pulsar kicks) for resonances occurring at p(r res ) = p c . 
In such a case, Eq. (15) reduces to V Vf ~ 0.7Vo (we use Y e ~ Y v ~ ©(lO -1 )). The mass of the sterile neutrino 
(m„ a rn Vf ) is derived through Eq. (12) for small mixing angle and p ~ 20MeV. One gets m Va ~ few keV, according 
to Ref. [22]. Eq. (38) then becomes 


E ~ 2 x 10 _ 2 a; res 77 eV. (40) 

Eqs. (24) and (40), and x res ~ 0(1) (as we shall determine below) yield 

|<5cr| < 4 x 10 _lo ? 7 , |<5ar| < 1.6 x 10 _ 2 ? 7 eV. (41) 

The constraints (5) provide rj < 10~'. 
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FIG. 1. The parameter fx, solution of Eq. (42), vs the parameter 77 . a is defined as a = r//10 4 , and, in the small panel it 
varies from 1 to 1 CP 3 . 

Eq. (28) admits the solutions x res = 1 and x res = \x/(\x — 1). The first solution is not compatible with (39) since 
it implies 77 = £ 2 A C /Ar + 3, so that no real solution there exists for e. The second one inserted into Eq. (39) yields 


£ = 


I A r (3 + ??) 
A c 


~ 0.82. 


Finally, the parameter /i is determined via Eq. (33) 


X 




l- L - 1 


— x(i) + 


2 A c 2Ar(3 77 ) 


= 0 , 


(42) 


from which one gets 77 ~ 17.21 (hence x res ~ 1.1 ~ 0(1)). Fig. 1 shows the values of the parameter /i for varying 77 . 
As one can see, fx approximates to 17.21 as 77 < 10~ 7 . 

From the above results, one infers 


T(r res ) ~ 1-2 T c , r res ~ 0.9 r c , i? s ~1.4r c . 


Let us now discuss the obtained results in relation to the velocity distribution of pulsars. At the moment, the 
statistical analysis of pulsar population neither support nor rule out any model/mechanism proposed to explain pulsar 
kicks. This is essentially due to the lacking of correlation between neutron star velocity and the other properties of 
neutron stars (see for example [35] for a general discussion and references therein). 

Concerning the coefficients for the Lorentz violation, it is convenient to adopt a standard inertial frame with respect 
to which experimental measurements of these coefficients are reported [36]. In this frame the coefficients a l and 
(hence 8 and 8c^) point along all possible orientation 2 . Usually, it is assumed that Lorentz violating parameters 
are independent on position [11 13]. Nevertheless, as discussed by Kostelecky, Lehnert, and Perry in [8], they may 
also exhibit a space-dependence, such that they may take arbitrary different values for pulsars at different places. 

These considerations entail two basic features: 

• A priori there is no favored direction for the Lorentz violating coefficients. This results in a uniform distribution 
of pulsars, in agreement with statistical analysis performed in Refs. [17,35,37]. 

• The mechanism proposed in this paper does not require any correlation between the velocity V and other 
physical parameters of pulsars, such as for example, the magnetic field B. The reason is essentially due to the 
fact that in the Standard Model Extension, the coefficients 0 ^ and Cj" describe the most general renormalizable 
effects that are possible in the gauge-invariant neutrino sector, no matter what are their origin at the Planck 


2 In particular, in the papers [ 11 , 12 ] it has been studied the general theory with a £ and which include all orientations, 
whereas in [13] the a)* coefficients point along the north direction, and care isotropic. 
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scale (they can be regarded as vacuum expectation values of tensor operators arising from the spontaneous 
breaking mechanism). As a consequence, a correlation among the coefficients and and the parameters 
characterizing the pulsar properties, in particular the magnetic field, is not necessarily expected. Such a result, 
indeed, seems to be corroborated by recent analysis of Refs. [17,35,37] in which it is pointed out the apparent 
lack of evidence in favor of the B — V correlation 3 . 

In conclusion, the origin of pulsar velocities has been studied in the framework of Standard Model Extension. The 
effective Hamiltonian for the time evolution of neutrinos (see Eq. (4)) gives rise to a term of the form ~ (511 • p. 
The relative orientation of neutrino momenta with respect to (511 leads to the fractional asymmetry necessary to 
the generate the observed pulsar motion. Future observations and statistical analysis on the velocity distribution of 
pulsars might provide a scenario in which the Lorentz and CPT invariance violation, as suggested by Kostelecky, 
might be tested on astrophysical scales. 
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APPENDIX: THE POLYTROPIC MODEL 


The aim of this Appendix is to recall the main topics of the polytropic model which describes the inner core of a 
protostar. To this end, we shall refer to the paper by Barkovich, D’Olivo, Montemayor and Zanella [23]. The relevant 
equations for the description of an isotropic neutrinosphere are the equation for the hydrodynamical equilibrium 


dP(r) , . GM(r ) 

—j — = -P{r) - 2 — > 

dr r 

where M(r) = 47r / Q r dr'r' 2 Pt{t') ( pr is the total density of matter), the equation for the energy transport 


(Al) 


and the equation for the flux conservation 


F(r) 


1 1 dT 2 

36 np(r) dr 


F(r) 


L c 

Anr 2 1 


(A2) 


(A3) 


where L c is the luminosity of the protostar. The equation of state with adiabatic index T is [34,23] 

P(r) = Kp v , K = - Tr 


m n p , 


1/3 


(A4) 


where T c and p c are the temperature and the matter density of the core. From Eqs. (A4) and (Al) one gets 


dp T 1 Ar r c Pc 1 Af (r) 

dr M c r 2 


(A5) 


3 It should be noted that the analysis of Refs. [17,35,37] does not imply that the Kusenko-Segre mechanism for generate the 
pulsar kicks, as well as those mechanisms that also rely on magnetic fields, do not work. As pointed out by Kusenko in [20], 
such mechanisms do not predict a B — V correlation since the relevant magnetic field for the kick velocity is the magnetic 
field inside the hot neutron star during the first seconds after the Supernova collapse. Astronomical observations allow to infer 
the surface magnetic field some millions of years later. These two fields are not trivially correlated each other because of the 
complex evolution of the magnetic field during the cooling process of neutron star, which lead to a final (surface) magnetic field 
whose configuration is different from the initial one, i.e. few seconds after the onset of the Supernova (see [20] for a detailed 
discussion). 
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where M c is the core mass. A function which represents an extremely good approximation for p is given by [23] 


P T 1 (0=dc 1 A r - l) m(r) + 


(A6) 


where m(r) = p + (l — p)r c /r. Notice that m(r c ) = 1. The parameter p is determined by using the equation p{R s ) =0 
( R s is the radius of the star), i.e. 


T = 


R s 


A r {R s - r c ) R s 


R s 


Rs - r n 


Eq. (A6) can be easily recast in the form 

P r - 1 ( r )=^ 1 


i'(-Y 

\ r / 


+ b'— + c' 
r 


where the coefficients a',b',c ' are defined in (29). Eqs. (A3) and (A2) allow to infer the temperature profile 

dT 2 3 4 5 6 7 8 9 nL r 


dr 7rr 2 ’ 

with p given by (A8). For T = 4/3 the solution of the integro-differential equation (A9) is 

T(r)=T c V2A c [x(r c /r)-x(l)]+l, 

with x( r c/i r ) defined in (34). 


(A7) 


(AS) 


(A9) 


(A10) 
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